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I. INTRODUCTION 

Diffractive scattering was imported into high energy strong interaction to interpret the t 
distribution of elastic hadron-hadron scattering at small |t| region |]l|. In the framework of the 
Regge theory, diffractive scattering is assumed to take place through the exchange of an object 
with the quantum numbers of the vacuum (the Pomeron) p[ , Although, the Pomeron plays a 
particular and very important role in soft processes in hadron-hadron collisions, the true nature 
of it and its interaction with hadrons remain a mystery. 

Due to the nature of color singlet exchange, the diffractive scattering events are characterized 
by a large rapidity gap, a region in rapidity devoiding of hadronic energy flow. Events containing 
rapidity gap and jets were first observed by UA8 Collaboration at CERN |^, opening the field 
of hard diffraction. The hard diffraction is also observed and studied by various experiments 
[^-^. Although it is still difficult to understand the soft diffractive scattering in the framework 
of QCD, much progresses in understanding the nature of hard diffraction are made in the 
light of great theoritical and experimental efforts |p. As an attempt of understanding hard 
diffraction on the parton level, Ingelman and Schlein [|lO| assumed that the Pomeron, similar 



to the nucleon, is composed of partons, mainly of gluons, and the hard diffractive processes 
can probe the partonic structure of the Pomeron. In fact, the Ingelman-Schlein model for 
hard diffraction is based on the assumption of hard diffractive scattering factorization together 
with Regge factorization for the Pomeron exchange. But, comparison the total pp{p) single 
diffraction cross section data []TT| with predictions based on the Regge factorization for the 



Pomeron exchange indicates the breakdown of the Regge factorization [0. Furthermore, as 
pointed by Ingelman in |P, there are conceptual and theoretical problems with the Ingelmal- 
Schlein model. First, the Pomeron is not a real state, it is only a virtual exchanged spacelike 
object, the concept of a structure function is then not well defined. In fact, the factorization of 
diffractive structure function into a Pomeron flux and a pomeron structure function cannot be 
uniquely deflned since only the product is an observable quantity. Second, the Pomeron-proton 
interaction is soft, its space-time scale is larger compared to that of the hard interaction, so the 
Pomeron can then not be considered as decoupled from the proton, and in particular, is not a 
separate part of the QCD evolution in the proton. 

For the lepton induced hard diffractive scattering processes, such as diffractive deep inelastic 
scattering (DDIS) and diffractive direct photoproduction of jets, there is indeed a factoriza- 
tion theorem which has been proven by Collins ||13[ recently. It is the same as factorization 
for inclusive hard processed, except that parton densities are replaced by diffractive parton 



densities |T^|T^. The diffractive parton densities (the primary non-perturbative quantities 
in hard diffraction) obey exactly the same DGLAP evolution equations as ordinary parton 
densities, and have been extracted from HERA data on DDIS and on diffractive photonproduc- 



tion of jets [0. This factorization theorem also establishes the universality of the diffractive 
parton distrubutions for those processes to which the theorem applies, hence justifles, from 
fundamental principles, the analysis of ZEUS and HI Collaboration p,^ on hard diffraction, 
provided the term " Pomeron " used in these analyses is as a label for a particular kind of 
parametrization for diffractive parton densities, and an indication of the vacuum quantum 
numbers to be exchanged. In the light of hard diffractive factorization proven by Collins [|13[] , 
the Ingelman-Schlein model for hard diffraction assumed further the factorization of diffractive 
parton densities into a universal Pomeron flux and a Pomeron parton densities, but such further 
factorization is unjustiflable. 

In this paper, we present a study of another hard diffractive process, J/ip + jet diffractive 
production in the direct photon process at HERA 



i{q) + p{pi) -^ p{p'i) + J I HP) + jet + X, (1) 

where q,Pi,p'i, and P are the momenta of incoming photon, proton, outgoing proton, and 
J/ip respectively. We base on the factorization theorem for lepton-induced hard diffractive 
scattering and the factorization formahsm of the nonrelativistic QCD (NRQCD) for quarkonia 
production. Using the diffractive gluon distribution function extracted from HERA data on 
DDIS and diffractive dijet photon production, we show that this process can give valuable 
insights in the color-octet mechanism for heavy quarkonia production. This process was briefly 



discussed in |[T^ where the relativistic corrections to photoproduction of J/ijj were studied in 
the so-called color-singlet model (CSM) |]T^. 

Our paper is organized as follows. We describe in detail our calculation scheme in Sec. II, in 
which, a brief introduction of the heavy quarkonium production mechanism, the hard subpro- 
cesses for J/ip + jet production and a summary of the kinematics related to the J/ip^s z and 
Pt distributiones are included. Our results and discussions are given in Sec. III. 



II. CALCULATING SCHEME 

Based on the factorization theorem for lepton-induced hard diffractive scattering and the 
factorization formalism of NRQCD for heavy quarkonia production, J/ip + jet diffractive pro- 
duction process is shown in Fig.l . This process consists two steps. First, the almost point-like 
cc pair with large Pt is produced in the hard subprocesses, then J/ip particle is produced from 
this point-like cc pair via soft interactions. 

A. Heavy quarkonium production mechnism 

Heavy quarkonium production in various processes has been the focus of much experimental 



and theoretical attention during the past few years ||I9[. This is mainly due to the observation 
of large disrepancies between experiment measurements of prompt and direct J/ip production 
and ip' production at large Pt at the Collider Detector Facility (CDF) at the Fermilab Tevatron 
and the calculations based on CMS pO[| . 

The first major conceptual advance in heavy quarkonium production was the realization that 
fragmentation dominates at sufficiently large Pt ||21|] which indicates that most charmonium at 



large Pt is produced by the fragmentation of individual large Pt partons. The fragmentation 
functions are calculated in CSM. Including this fragmentation mechanism indeed brings the 
theoretical predictions for prompt J/ip production at the Tevatron to within a factor of 3 of 
the data [^. But the prediction for the ip' production cross section remains a factor of 30 



below the data even after including the fragmentation contribution (the ^/''''surplus" problem). 
Furthermore, the presence of the logarithmic infrared divergences in the production cross sc- 
tions for P-wave charmonium states and the annihilation rate for Xcj ~^ QQg indicate CSM 
is incomplete. All these indicate that important production mechanism beyond CSM needs 
to be included ||2^ p4|. So the color-octet mechanism ||19[ is proposed which is based on the 



factorization formalism of NRQCD |2^ ||26|. Contrary to the basic assumption of CSM, the 
heavy quark pair in a color-octet state can make a transition into physical quarkonium state 
through soft color interactions. 

Although the color-octet fragmentation picture of heavy quarkonium production p4| has 



provided valuable insight, the approximation that enter into fragmentation computations break 



down when a quarkonium's energy becomes comparable to its mass. The fragmentation pre- 
dictions for heavy quarkonium production are therefore unreahable at low P^. Based upon the 
above several recently developed ideas in heavy quarkonium physics, Cho and Leibovich p7 



identify a large class of color-octet diagrams that mediate quarkonia production at all energies, 
which reduce to the dominant set of gluon fragmentation graphs in the high Pt limit. By 
fitting the data on heavy quarkonia production at the Tevatron, numerical values for the long 
distance matrix elements are extracted P7| , p5| , which are generally consistent with NRQCD 



power scaling rules [2S 



In order to convincingly establish the color-octet machanism, it is important to test whether 
the same matrix elements be able to explain heavy quarkonium production in other high energy 
processes, such as inclusive J/ip production in e~^e~ annihilation and inelastic J/ip photopro- 
duction at HERA et al. . Other machanism for heavy quarkonium production can be found in 



30,31 



B. The hard subprocesses for J/ip + jet production 

J/ijj is described within the NRQCD framework in terms of Fock state decompositions as 

I J/^) = 0(1) \ccfS['^]) + Oiv)\ccfPi%) 

+Oiv')\cc['SJ,%) + Oiv')\cc['si'%g) 

+Oiv')\cc['p!j'%g) + ---, (2) 

where the cc pairs are indicated within the square brackets in spectroscopic notation. The 
pairs' color states are indicate by singlet (1) or octet (8) superscripts. The color octet cc state 
can make a transition into a physical J/ip state by soft chromoelectric dipole (El) transition(s) 
or chromomagnitic dipole Ml transition(s) 

{ccr^'Lf]^J/^. (3) 

NRQCD factorization scheme ||2^ has been established to systematically separate high and 
low energy scale interactions. It is based upon a bouble power series expansion in the strong 
interaction fine structure constant a^ and the small velocity parameter v. The production of a 
{cc)\^^~^^Lj ' ] pair with separation less than or of order ^ can be calculated perturbatively. 
The long distance effects for the produced almost point-like cc to form the bound state are 
isolated into nonperturbative matrix elements. Furthermore, NRQCD power counting rules 
can be exploited to determine the dominant contributions to various quarkonium processes [ p9[| . 
For direct J/ip production, the color-octet matrix elements, (0|(9g pS'i]|0), (OlCg ["^5'o]|0) and 
(0|Og pPj]|0)/m^ are all scaling as mlvl. So these color-octet contributions to J/ijj production 
must be included for consistency. 

On the partonic level, J/ip + jet production are composed of the photon gluon fusion, which 
are sketched in Fig. 2. These are 

^ + g-.g + {cc)['s['\' S['^] , Fig.2.(a) - (f); 
l + g^g+{cc)['SJ,'\'P!f^] , Fig.2.(a) - (h). 

The quark initiated subprocesses(7g(g) channel) are strongly suppressed and will be neglected 
further. 

The color SU(3) coefficients are given by 



(3«; 3j|l) = 5j,/V3, (3^; 3j|8a) = v^T/,. 



(4) 



Some identities involving the traces of the color matrices are useful when the matrix elements 
are squared and spin, color averaged: 

iTi.(T'r=)Ti.-(r'r)-L-L = ^. 
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Where 



p(a)-(c) _ y/2Tr (T'^T^T^) 
p(d)-(f) _ y2Tr(T"T"T^), 
F(g)-(h) ^ y2Tr(T'^T'^)/rffee 

„ Jabc- 



(6) 



.^Arlc(8) 3 p(8)i 



Note that the hard subprocesses that producing color-octet {cc)[^Sq , P) ] pair involve 
triple-gluon coulping (Two of them are external gluons). This situation implies that we cannot 
use the trick of replacing the sum over gluon polarization states by 



^e„(A)e;,(A) ^ -go 



(7) 



when we use the Feynman gauge for the gluon propagator |^. One way to obtain the cor- 
rect result is insist that the polarization states of the two external gluons are physical {i.e., 
transverse). This is accomplished by the projection 

A 

Qaa' 



(nk) 



+ 



{nky 



where n is an arbitrary 4-vector and k is the gluon 4-momentum. This is analogous to using 
an axial (z.e. physical) gauge for the gluon propagator. Our choice for the 4-vector n is n = g2 



for the incoming gluon and n = gi for the outgoing gluon, gi,g2 are the 4- momenta of the 
incoming and outgoing gluons respectively. 

With all ingredients set as above, we obtain the spin, color average- squared matrix elements 
E|M(7 + g ^ {cc)[^^~^^Lj ' ] + g)\'^. Our analytic expressions are consistent with [^] and the 



results calculated from Eqs. (11),(15) and appendix B of |Q, but different from the expressions 
of [0. For example, the expression of E|M(7 + g ^ (cc)[^S'o ] + g)\'^ of [Q is different from 
our analytic expression by a term 



6su{eecg 
i{s + u) 



2n2 



(9) 



which is exactly the contribution of photoh-ghost scattering diagram [i.e. the contribution of 
the unphysical polarization states of the external gluons) and should be subtracted from Eq. 
(Al) of ig. Here 



s = {q + 9i)\ i={q-Pf 



u 



(91 - Pf- 



(10) 
+ g)\'^ for J/ip pro- 



In appendix, we give explicit expressions of S|M(7 + g 
duct ion. 

Known the spin, color average-squared matrix elements of the hard subprocesses, in the 
NRQCD framework, we obtain the subcross sections. The color-singlet photon-gluon fusion 
contribution to J/-?/' -|- jet production is well known ||3^ : 
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where (0|(9^ ['^S'i]|0) is the color-singlet matrix element which is related to the lepton decay 
width of J/ip. The average-squared matrix element S|M(7 + g ^ (cc)pS'j ] + g)\'^ can be 
obtained from S|M(7 + g ^ {cc)^Si ] + g)\'^ by taking into account of different color factor. 
The color-octet (cc)[^S'j ] contribution is 



da 
di 



b + 9 



15 



IGtts^ 6 
1 , 
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2Amr 



S|M(7 + (7- 



J/i^ + g] 
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:i2) 
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rs p(8)i 



The color-octet {cc)[^Sq ] and {cc)['^Pj ] contributions are 

da , 



dt 



[^ + g^{cc)['si,'^]+g^J/^+jet] 



lc(8)l 



IGtts^ 



E\Mi^ + g^icc)['Sn+9)\ 



'A'So 



^ [7 + ^ -- (cc) fP^f^] +g^J/^ + jet] 
at 



16iiP J 
X ^{0\Oi/^['Pom, (13) 

where the heavy quark spin symmetry 

(O|08^/'^pP^]|O) = {2J+l){0\Oi^^['Po]\0) (14) 

is exploited. 

C. The z and Pt distributions of J/ip 

Now we consider the z and Pt distribution of J/ip produced in process Eq.(0). Based 
on the factorization theorem for the lepton induced hard diffractive hard scattering |T^, the 



differential cross section can be expressed in terms of a diffractive parton distribution [p^pTS 
as 

dfg}"{xi,XTP,t,fi) 

da = — ^^^ — ; ; dxjpdt ■ 

dx-pdt 

^^^^l±^^M±^dx,di. (15) 

dt 

Here 

<"':"7-'-"'d.. (16) 

dxjpdt 

represents the probability of finding in the proton a gluon carrying momentum faction Xi, while 
leaving the proton intact except for been diffractively scattered with the momentum transfer 
(xp,t). Xp is the fractional momentum loss of the diffracted proton, i.e., Xp ~ {piz ~ Piz)/Pz, 
and t is the invariant momentum transfer for the diffracted proton, i.e., t = {pi — p'l)"^- 
We now consider the kinematics. It is convenient to introduce the variable 

..^^. (17) 

Pi -q 

In the proton rest frame z = E^/ E^. In terms of z and Pt (the transverse momentum of J /ip), 
the Mandelstam variables and xi can be expressed as the following: 
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Here, s^p = (g + p^- The double differential cross section is 

da /-^ipmax fO dfyj{xi,xjp,t,^) 

dxjp I dt — ^^-^ 



dzdPr J XI J-i dxTpdt 

daj'y + g ^ J/ip + jet) xii 
di zPt 



where the Jacobian can obtain from Eq.(|T8| 



A^) = '"^"t ■ (20) 

^zPt' z{1-z)s^p ^ ' 



The allowed regions oi z,Pt are given by 



< Pr < \/(l - z){zXTPraa^Syp - m^), 



with 



■^•^Pmax'S'yp 



■^IPmax'^^p + ^^ip 






■^•''Pmax''7p 



'7P ^ ""i/j 



In order to suppress the Reggon contributions, we set a;pmax = 0.05 as usual. 

III. NUMERICAL RESULTS AND DISCUSSIONS 

For numberical predictions, we use mc = 1.5GeV,A4 = 235MeV, and set the factor- 
iztion scale and the renormalization scale both equal to the transverse mass of J/ip, ^-e., 
^^ = m^ = {ml + P|) . For the color-octet matrix elements {0\Oi^^[^Si\\0), (OlOi^'^l^SoWO) 

and (O|08 pPo]|0) we use the values determined by Beneke and Kramer |2^ from fitting the 
direct J/ip production data at the Tevatron using GRV LO parton distribution functions PB 



{0\Oi^^[^Si]\0) = 1.12 X lO'^GeV' 



^lj|w/ - 

3.5 

m^ 
= 3.90 X IQ-^GeV^. (23) 

Since the matrix elements (OjCg [-'^S'o]|0) and (OjCg pPo]|0) are not determined separately, 
we choose 

{0\Oi^^[^So]\0) = 1.0 X lO-^GeV^, 
{0\Oi^^[^Po]\0)/ml = 8.3 X 10-3Gev^ (24) 



allowed by Eg. (]23|) , and in accordance with NRQCD power counting rule. The value of the 
color-singlet matrix element is taken to be (0|(9i p^iJlO) = 1.16GeV^. The diffractive gluon 
distribution function 



dx-pdt 



(25) 



at /i^ = 4GeV^ has been extracted from data on DDIS and on diffractive photoproduction of 



jets at HERA in [16|, 



dfg'//{xi,xjp,t,fi'^ 



4GeV^ 
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-2a{t) 
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ag{l -xi/x-p), 



(26) 



Tl. Here po = l.SGeV" 



4.5 ± 0.5. We use the central value of ag for numberical calculation. 



and is evolved in /x^ according to the DGLAP evolution equations 

a{t) = 1.14 + 0.25t, Oj 

In order to suppress the elastic J/tp photoproduction (with or without proton dissociation) 

contribution and higher-twist corrections, a Pt cut Pt > l.OGeV is imposed. 

In Fig. 3 we show the z distribution ^ of J/ip produced in process Eq.(|l]) at HERA at a 
typical energy ,/s^ = lOOGeV and with a Pt cut Pt > l.OGeV. The thin dash-dotted and short 



dashed lines represent the color-singlet cotribution and the sum of the color-singlet and color- 
octet contributions respectivly. We observe that for z > 0.3 the color-singlet contribution is 
about order of 10°nb, so this process can be studied at HERA with present integrated luminosity. 
For completeness, in Fig.3, we also show the z distribution of J/ip inelastic production 'j + p-^ 
J/ip + jet + X in the same kinematic region (upper part), in this case, the thick solid and 
dotted lines curves are the sum of the color-singlet and color-octet contributions and the color- 
singlet contribution respectively. In both case, including the color-octet contribution, the z 
distribution is strongly enhanced in high z region due to the gluon propagator in the color- 
octet channel (Fig.2.(g)-(h)). In inelastic case, this behavior of rapid growing at high z does 



not agree with the HERA data []37| , p8 . 
In Fig. 4, we show the Pt distribution 



at HERA with 



lOOGeV and a z cut 



da 

z < 0.8. The code for the curves are the same as Fig.3. We find that for J/ip + jet diffractive 
photoproduction in process Eq. (|l|) the color-octet contributions are larger than the color-singlet 
contribution for Pt > 2.6GeV. 

In Fig. 5, we show the J/ip + jet diffractive photoproduction cross section as a function of 
■,ys^ with the cut, Pt > l.OGeV, z < 0.8. We observe that in the energy region considered the 
color-singlet contribution is about order of 10°nb, so this process can be studied at HERA with 
present integrated luminosity, and can give valuable insights in the color-octet mechanism for 
heavy quarkonia production. 

Recently, it was pointed out that the cross section with an additional z cut, say, z < 0.8, 
cannot be reliably predicted in NRQCD p4| , p9| , because the NRQCD expansion is singular 
at z = 1, only an average cross section over a sufficiently large region close to 2; = 1 can be 
predicted. The z distribution itself requires additional non-perturbative information in the 
form of so-called shape functions. These shape functions are also required to predict the Pt 
distributions with an additional upper z cut, say, z < 0.8, but not if z is integrated up to 
its kinematic maximum. So in Fig. 6, we show the the Pt distribution -^ at HERA with 



Js^p = lOOGeV and with z is integrated from its lower cut 0.2 to its kinematic maximum. The 
code for the curves are the same as Fig. 4. We find that for J/ip+jet diffractive photoproduction 
in process Eq.(|l|), without the upper z cut, the color-octet contributions are dominated in the 
whole Pt region considered, which exceed the color-singlet contribution by almost an order 
of magnitude. In Fig. 7, the J/ip + jet diffractive photoproduction cross section as a function 
of ^JSryp with z is integrated from its lower cut 0.2 to its kinematic maximum is shown. We 
observe that in the energy region considered the color-singlet contribution is in the region 
1.4nb < a'^'^^^' [singlet] < 2.4nb, while including the color-octet contributions, the cross section 
increased by about an order of magnitude. 

Experimentally, the nondiffractive and elastic (with or without proton dissociation) back- 
ground to the diffractive J/ip photoproduction we studied must be dropped out in order to 
obtain useful information about the heavy qarkonia production mechanism, this can be at- 
tained by performing the rapid gap analysis together with a large Pt cut, say, Pt > l.OGeV. 
If statisitics is not a limitation, it might be preferable to use cut Pt > 2GeV. 

In conclusion, in this paper. We present a study of J/ip + jet diffractive production in the 
direct photon process at HERA based on the factorization theorem for lepton-induced hard 
diffractive scattering and the factorization formalism of the nonrelativistic QCD (NRQCD) for 
quarkonia production. Using the diffractive gluon distribution function extracted from HERA 
data on diffractive deep inelastic scattering and diffractive dijet photon production, we show 
that this process can be studied at HERA with present integrated luminosity, and can give 
valuable insights in the color-octet mechanism for heavy quarkonia production. 
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APPENDIX: THE SPIN, COLOR AVERAGE-SQUARED MATRIX ELEMENTS 
In this appendix, we give explicit expressions for the spin, color average-squared matrix 
elements of photon-gluon fusion subprocesses. The results are obtained using the symbolic 
manipulations with the aid of REDUCE package. 
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Fig.l. Sketch diagram for J/f/i + jet diffractive production. 
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Fig. 2. Feynman diagrams for 7 + 5 — > cc[^^"'"-'ij ] + g subprocess. (a)-(f) 
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Fig. 3. z distribution da jdz for J/%lj+jet production at HERA at JSyp = lOOGeV, integrated 
over the J/ip transverse momentum Pt with a lower P^ cut P^ > l.OGeV. The thin dash-dotted 
and short dashed hues represent the color-singlet cotribution and the sum of the color-singlet 
and color-octet contributions to J/ip+jet diffractive photoproduction 7+p -^ J/ip+jet+p+X 
respectivly. The thick solid and dotted lines curves are the sum of the color-singlet and color- 
octet contributions and the color-singlet contribution to inelastic J/ip + jet photoproduction 
7 + p — > J/ip + jet + X respectively. 
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Fig. 4. Transverse momentum (Pt) distribution da/dPr for J/ip+jet production at HERA 
at y/s^ = lOOGeV, integrated over z with a upper z cut z < 0.8. The thin dash-dotted and 
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7 + p -^ J/ip + jet + X respectively. 



10' p-TT 



J2 



10 




10' - 



10" — 



TOT-inelastic 
CS-inelastic 
TOT-diffractive 
CS-diffractive 



I I I I I I 



_L 



100 



150 

s 



200 

(GeV) 



250 
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